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Abstract 

The equation of state for a two-dimensional hard-sphere gas is difficult to calculate by usual methods. 
In this paper we develop an approach for calculating the equation of state of hard-sphere gases, both for 
two- and three-dimensional cases. By regarding a hard-sphere gas as an ideal gas confined in a container 
with a multi-core (excluded sphere) boundary, we treat the hard-sphere interaction in an interacting gas as 
the boundary effect on an ideal quantum gas; this enables us to treat an interacting gas as an ideal one. We 
calculate the equation of state for a three-dimensional hard-sphere gas with spin j, and compare it with the 
results obtained by other methods. By this approach the equation of state for a two-dimensional hard-sphere 
gas can be calculated directly. 

PACS numbers: 05.30.-d, 05.70.Ce 

There are various methods for obtaining the equation of state of three-dimensional hard-sphere gases, e.g., 
the method of pseudopotentials 1 , cluster expansions [2 El- However, it is difficult to calculate the equation 
of state of two-dimensional hard-sphere gases. In this paper, for dealing with the two-dimensional hard-sphere 
gases, we develop an approach for calculating the equation of state of hard-sphere gases, both for two- and 
three- dimensional cases. This paper falls into three sections. In section 1 we first introduce the basic idea of 
this approach. For illustrating the approach and demonstrating its validity, we use this approach to calculate 
the equation of state for three-dimensional hard-sphere gases and compare the result with the standard result 
given by Lee and Yang [21 El- In section 2 we calculate the equation of state of two-dimensional hard-sphere 
gases. The paper is concluded in section 3. 

1 Hard-sphere gases as ideal gases with multi-core boundaries: An 
approach to two- and three-dimensional interacting gases 

By regarding a hard-sphere gas as an ideal gas confined in a container with a multi-core boundary, i.e., a 
container filled with small excluded spheres (Fig. 1), we calculate the grand potential for two- and three- 
dimensional interacting gases using the method developed for calculating the boundary effect on an ideal gas. 
The value of studying hard-sphere gases is: (1) Often a complete knowledge of the detailed interaction potential 
is not necessary for a satisfactory description of the system, because a particle that is spread out in space sees 
only an averaged effect of the potential jT|. (2) Experimental progress has been made in recent years in Bose- 
Einstein condensation in dilute atomic gases the hard-sphere interaction model is often used as an efficient 
tool for studying Bose-Einstein condensation, both for dilute atomic gases and liquid helium 5 , after Lee and 
Yang studied pioneered the study of interacting bosons 1^1- (3) The hard-sphere gas, as a simplified model, is 
of great value for investigating the more general theory of interacting gases and can be extended to some more 
general cases. As a valuable idealized model, the hard-sphere gas has been studied by various theories, such as 
the method of pseudopotentials 1. , cluster expansions |21 El , and second quantization \7\ . 

The hard-sphere gas is a simplified model of the interacting gas, which replaces the interparticle interaction 
by the boundary condition the wavefunction ?/; = on the boundary; such a boundary in the 3N-dimensional 
configuration space is equivalent to a tree- like hypersurface UEl. In a quantum hard-sphere gas, there are two 
interplayed effects: the effect of the statistics and the effect of the hard-sphere interaction. The interplay of the 
effects of statistics and the interparticle interaction often causes difficulties, e.g., it makes the quantum cluster 
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tEmaihdai wusheng@tju.edu.cn 
-fEmaihxiemiQtju. edu.cn 



1 



expansion considerably involved [H]: to find the Z-th virial coefficient for I > 2, one has to solve a quantum 
many-body problem. There are various approaches to hard-sphere gases. The scheme proposed by Lee and 
Yang, the binary collision method, is to separate out the effect of the statistics |3 El ; it allows one first takes 
care of the statistical aspect of the problem and then tackles the dynamical aspect of it , and the l-th virial 
coefficient can be obtained by solving two-body problems. In the method of pseudopotentials the boundary 
condition is approximately replaced by a non-Hermitian pseudopotential mH], i-e., the boundary condition is 
converted to a "potential" approximately. 

In this paper, by replacing a hard-sphere gas by an ideal quantum gas with a multi-core boundary (Fig. 1), 
and using the method developed for calculating the boundary effects, we converted the problem of a quantum 
hard-sphere gas (Fig. la) into a problem of an ideal quantum gas confined in a container filled with small 
excluded spheres (cores) randomly distributed (Fig. lb). In this approach, roughly speaking, the exchange 
effect and the effect of classical interaction (in this case, the hard-sphere interaction) are treated separately: the 
exchange effect on the imperfect gas is embodied in the exchange effect on the ideal gas; the effect of hard-sphere 
interaction is embodied in the boundary effect. The number of the cores is chosen to equal the total number of 
particles; the radius of the cores is chosen to equal the diameter of the hard sphere a (the scattering length of 
the hard sphere) since the area of a core felt by a gas molecule during the scattering process is Aira^. 




Fig. 1. Regarding a hard-sphere gas as an ideal quantum gas in a multi-core container, (a) A hard-sphere gas in which 
the diameter of the gas molecule is a. (b) An ideal quantum gas confined in a container with a multi-core boundary. 
The diameter of the core is 2a. 

Then, the task that to calculate an interacting gas is converted to a task that to calculate the boundary 
effect on an ideal quantum gas. In preceding papers we developed some methods for calculating boundary 
effects ^2 El- one which is based on the mathematical work given by Kac "W , is for two-dimensional 
cases; the other ^21 is for three-dimensional cases. However, to apply these methods to the present problem, in 
which the boundary is more complex, we need to extend them to more general cases |14| : the grand potentials 
for ideal quantum gases in three- and two- dimensional confined space can be expressed as 

^^- = 97m h5/2{z)-g-^ -2 h2{z) + g-^ /13/2 (^) , (3D) (1) 



h2{z)-g- 



In S ^ 4 - ''^ hAz)-g Y'- ^-'r ^3/. (^) , (2D) (2) 



where d denotes the dimension; A = h/\/^jrmkT is the thermal wavelength, g the number of internal degrees of 
freedom and for spins g = 2j + 1, and the function hcr{z) = z-^e^^i equals the Bose-Einstein integral 

ga{z) or the Fermi-Dirac integral fa{z) in boson or fermion case respectively. ^° and /i^ are the valuations 
of the containers 15 (/i° for outer boundaries and /i^ for inner boundaries): in three dimensions for a sphere 
Hz is the volume, /i2 half the surface area, and /^i twice the diameter; in two dimensions, for a disk /i2 is the 
area, and half the perimeter. However, this result is only valid for spinless distinguishable particles. To 
take the effects of spin and indistinguishability into account, we need to improve this result one step further. 
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The hard-sphere interaction is embodied in the boundary terms; the cross-section for colhsions between gas 
molecules is regarded as the surface area of the boundary cores. For distinguishable hard-sphere particles, the 
cross-section is Ana^; however, for indistinguishable ones, the cross-section is Swa^, i.e., the contribution from 
the indistinguishable hard-sphere interaction is twice the contribution from distinguishable ones. Moreover, for 
the case of particles with spin j, if we include only the s-wave contribution, there are weights (j -I- 1) / (2j -I- 1) 
for Bosons and j/ {2j + 1) for Fermions respectively for the unpolarized scattering of identical particles 
Thus the grand potential of an ideal quantum gas confined in the container illustrated in Fig. lb can be written 
in the following form: 



„ V uj 
ln^ = .9— /i5/2 W+.92^-— Y 



1 iiV47ra2 

2 A2 



h2 [z) 



-2N2a^ 
A 



^3/2 



(^) 



(3) 



where to = j + l for bosons and uj — j for fermions, a is the diameter of the hard sphere and N is the total number 
of particles in the system. Note that the boundary, of course, includes all cores. Strictly speaking, see Fig. 1, 
the volume of the system is the volume of the cubical container V minus the total volume of all cores; however, 
in view of the fact that for a dilute system the former is much larger than the latter and the contribution from 
the volume of the cores is only of the order [a/Xf" ^ we approximately use V as the volume of the system. In 
the second term the factor (1/2) NAtto^ is half the total surface area of all cores — the second valuation, which 
is much larger than the magnitude of y^/-^ in a macroscopic system; in the third term the factor 2 x 2a is twice 
the diameter of a core so 2N2a is the third valuation of the boundary (notice that 2N2a ^ V^^^). The factors 
2 and uj/ {2j + 1) reflect the effects of indistinguishability and spin. From Eq. Q, the relation N — zdlnE/dz 
gives a differential equation in N. Such an equation for N is solvable (see Appendix); however, its solution is 
in a complicated form. In the following, we present a series solution, which is easy to be compared with the 
results given in the literature, and, of course, it can be checked directly that such a series solution is consistent 
with the solution obtained by directly solving the differential equation for N. For seeking a series solution for 
A'^, we express A as a combination of the Bose-Einstein or the Fermi-Dirac integrals, ha{z), 



N = g 



V 



i<-n<\ 



«<'?<A 



e<';<A 



(4) 



where the sum is taken by all possible h„{z) and the subscripts of the coefficients are taken in ascending order 
for avoiding repetition, and substituting Eq. Q into Eq. accurate to (a/ A) , we obtain 



, _ y I , UJ 



X '2j + l 



A^h3/2hf_, + ^ B^^hs/2h^h„ + ^ Q,,A^3/2^5S'iA 



C<,,<A 



(a) '^2jTT 



?<';<A 



S.<ri<X 



(5) 



At v<a 

Deriving N by using the relations N = zdln'E./dz and zdh^/dz = /io-_i, after ignoring higher-order terms, we 
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have 



N = g 



a LJ 



(a) 2-^27TT 



u<,a 



1><<T 



hi/2hi,ha + ^ B'^,^h^i2hi,^iha + ^ Bl„h^i2h^h„^i 



and comparing the expression of with Eq. Q), we obtain the coefficients: ^43/2 = 1, B'^/2 3/2 ^ ^4^1^/ (2j + 1), 

^13/2 ~ ^1/2 2 ~ — STrgcj/ (2j + 1), — 2C"^^2 3/2 3/2 ~ ^^9^ ^1 C^-j + -'^)]^' other coefficients 

equal 0. Then, we achieve the equation of state, Eq. © 



^3/2 ^ ^ (2i + 1) 27r/l3/2/l2 



(6) 



Here we have substituted the relation g = 2j + 1 and replaced the Bose-Einstein integral (z) by her (z), which 
includes the Fermi case. This is the result for an ideal gas confined in the container. Fig. fb, and this is also 
the result what we want to obtain for a hard-sphere gas. Fig. la. 

Alternatively, we can also obtain this result in a more rigorous way. For seeking the series solution for N, 
we can also express N as 



00 00 2 



(7) 



In the following we take the Bose-Einstein case as an example, and the result for the Fermi-Dirac case can 
be achieved by a very similar procedure. Substituting Eq. Q into Eq. ^ and expanding the Bose-Einstein 
integral as a series by use of g^ = J2^i /I" , we have 



A3 1 ^ Z5/2 A ^27 + 1 
1=1 ■' 

-I 00 . 

J + 1 A 
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E 
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The relation N =^ zdlnR/dz, after ignoring higher-order terms, gives 

I °° I -11 °° /I 

I . 2 „ 7 + 1 a ^ . A 

fi,i=i 
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Comparing with Eq. , we obtain 
1 



4 — n 
~ .,3/2 ' 
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E 
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where the summation J2^a should be regarded as when b < a. Substituting these parameters into Eq. © gives 



InS 



1=1 



J + 1 
2j + l 



a \ ^ 



1 1 



A ^ ^3/2 /(3/2 



(ay ( j + 1 ^ 2 I L 

\XJ \ '^2j + l ^ ;3/2 1/2)^3/2' 



^ 1 1 

^ 1^ 12 3/2- 



.i-i+i 



Performing the summations and substituting the relation g = 2j + \, we again achieve the equation of state, 
Eq. ®. 

Our result Eq. (jSJ compares well with the result based on the binary collision method given by Lee and 
Yang El, 



(9) 



where F (z) = YlTs t=i irst)~^^'^ (r + s)^^ (r + t)~^ z''+''+*, "+" for bosons and "— " for fermions. The first- 
order contributions, a/ A, are the same; the difference only appears in the second term of the second-order, 
(a/A) , contribution. To show this difference clearly, as an example, we compare them in the low-temperature 
and high-density Fermi case, which is the most interesting case — the quantum case. The asymptotic forms for 
large z are 8F {—z) ~ —2.10 (In z) and —2nh^/2h2 ^ —2.36 (In z) ' . For Bose cases, at the high temperatures 
and low temperatures, our result is consistent with Lee and Yang's. At low temperatures and high densities, 
it is physically meaningless to compare the results Eqs. © and © since they are both invalid due to the 
existence of the critical point; in fact in this case the perturbation theory in a makes no sense j^l I17|. Notice 
that the methods given in Ref. |3| and in this paper are both approximate methods giving fugacity series of 
grand potentials. 



2 The equation of state for two-dimensional hard sphere gases 

The two-dimensional hard-sphere quantum gases which are difficult to calculate by usual methods can also be 
calculated by this approach. Like the three-dimensional case, the grand potential for the two-dimensional cases 
can be expressed as 



InS = g 



1 



S-2- 



Nna'^ ] h2- g 



2j + l 

g^h2 - g2 Nh'./o ~ q2 

^ ^ 2j + l\2 ^'^ ^ 2j + l 



^ 2^^N-27Ta]h 
2 A V 2j + 1 2 



3/2 



-j nNh2. 



(10) 



It should be emphasized that, in the two-dimensional case we take the area of the cores into account since it pro- 
vides a second-order, (a/A)^, contribution (Recall that in the three-dimensional case we ignore the contribution 
coming from the volume of the cores since the volume contribution is of the order (a/ A) ). 
By analyzing the grand potential Eq. H10|l . we express the particle number TV as 



N = g^ 



e<r,<A 



1><<T 



C<,,<A 



?<';<A 



(11) 
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Substituting N into Eq. Hl()|l . we have 
„ S S bj a 



C<r,<A 



S* Ufa 



A2 "2j + l VA 



S* w /a 



«<'J<A 



Af_,h2hfj^ + y] Bi,„h2Kh„ + C^r,xh2h^hrjh), 



(12) 



By = Z(91n^/9z, after Ignoring higher-order terms, we have 
N^g—hi 



S 
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Comparing the coefficients with the corresponding coefficients in {HJ, we have 



■2j + l' ^'^/^ "2j 

A" = O.B" -2719-^,5" = -2715-^, 
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= 3 7r.g 



C-1,3/2,3/2 - \^9i^-^ 
Then we obtain the equation of state for two-dimensional gases (17 = 2j + 1): 

- (2j + l)/j2+c^(2j + l) 



--7r/ii/i3/2 -t- j (^cjvr /io/i3/2 + t^TT hii2hih^/2 - 27r/ii/i2 
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Alternatively we can also obtain the equation of state in a more rigorous way. For seeking a series solution 
for TV, we expand N as 



Substituting into Eq. pU|l . taking the Bose case as an example, we have (a; = j + 1) 
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Then we can calculate N: 
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Comparing with (|14|l we obtain the coefficients: 
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J + 1 
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EE 
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Substituting these coefficients into Eq. H15() . we have 



2^5 



J + 1 



2j + l ttM(fc-M) 



(16) 



InS 



= .9^E 



S" a 
A 



/ + 1 j ^3/2 



^+1 

2j 



2 / j + l 

TTg- 
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2 oo 



E - 
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Therefore, following a similar procedure to that which led to Eq. JB)), we obtain the equation of state for the 
two-dimensional hard-sphere gas, Eq. again. 



3 Conclusions 

In conclusion, by using a method developed for calculating the boundary effect we treat a hard-sphere gas as 
an ideal quantum gas in a multi-core container, and calculate the grand potentials for hard-sphere gases in 
two and three dimensions. The result shows that by only replacing the moving gas molecules with a fixed 
boundary, we can obtain the results almost accurate to the second order, (a/A)^. This is because in this 
approach we use the s-wave scattering cross-section as the area of the cores. When the p-wave contribution 
is taken into account, one should consider the effect of the motion of the cores. Such a simplified model for 
interacting gases allows us to calculate the hard-sphere gas in a simple manner, and can be used to calculate the 
boundary effects in interacting quantum gases since in this treatment the hard-sphere interaction is regarded 
as a boundary effect, which we will discuss elsewhere. One advantage of this approach is that it can be used to 
calculate the hard-sphere gas confined in a container with the Dirichlet boundary condition since in this paper 
the method used to calculate the boundary effect is valid both for periodic and Dirichlet boundary conditions; 
while, e.g., in the pseudopotential method the periodic boundary play an essential role in the calculation since 
with such a boundary condition the momentum is conserved jH|. For seeking more higher-order solutions, we 
need to take into account other corrections, e.g., the correction caused by the movement of the cores, etc. 
In such a consideration the grand potential cannot be expanded in powers of a/A as the higher-order result 
obtained by using the method of pseudopotentials JHI- This method can be extended to more general cases: 
by replacing the diameter of hard-spheres with the scattering length, we can obtain a generalized result which 
can be applied to other kinds of interacting gases. Moreover, this result still can be carried one step further. 
The expansion coefficients can be explained as form factors. In the case of hard-sphere interactions, such form 
factors are constants. If we introduce a set of parameter-dependent form factors, the method can be generalized 
to describe more complex interacting gases. Such a treatment, which will be discussed in detail elsewhere, is 
equivalent to replacing the hard sphere by a soft sphere which has inner structures. Moreover, this approach 
can also be directly applied to quantum gases or liquids in porous media. Porous media can be regarded as 
containers with many holes (cores); in this case the particle number N in Eq. Q should be replaced by the 
number of such holes (cores) which are actually real boundaries of the container. 



7 



We are very indebted to Dr. G. Zeitrauman for his encouragement. This work is supported in part by the 
science fund of Tianjin University and Education Council of Tianjin, P. R. China, under Project No. 20040506. 
The computation of this project was performed on the HP-SC45 Sigma-X parallel computer of ITP and ICTS, 
CAS. 



4 Appendix [il9j| 

As mentioned above, the equation for N is solvable; in this Appendix we give the solution for this equation. 
From Eq. the relation N = zdlnE/dz gives a differential equation for N, 



dz uz {^h3/2 + 
The solution of this equation is 



N- 



N 



J ~p(z)dz 



where 



(2j + 1) y,h^/2 

^ UJZ (^/l3/2 + ^h2) 

Ip(^)d^q{z)dz 



= 0. 



(17) 



(18) 



p{z) = 



L^z{^hy2 + ^h2) ' 
(2j + 1) ^hy2 



3/2 



However the form of this solution is complex, and it is difficult to covert this solution to a simple and more 
useful form. Of course, the solution given by Eq. (|18() is consistent with the series solution given in the paper. 
To check this, we can, for example, expand Eq. (|18|) and compare it with the expansion of the result given in 
the paper. When z < 1, the expansion of Eq. H18|) and the result of N in the paper both are 



iV=(2j + l)^z+(2j + l)^ 



V2 , a ^ 

4a; 47rcLi 

4 A 



z^ + 



They are the same. In a word, what we do in the paper is to find a series solution for the differential equation 
of iV. 



References 

[1] HUANG K., YANG C.N., Phys. Rev. 105 (1957) 767; 

HUANG K., YANG C.N., and LUTTINGER J.M., Phys. Rev. 105 (1957) 776. 

[2] LEE T.D. and YANG C.N., Phys. Rev. 113 (1959) 1165. 

[3] LEE T.D. and YANG C.N., Phys. Rev. 116 (1959) 25. Before the publication of this paper, the results had 
been reported in LEE T.D. and YANG C.N., Phys. Rev. 105 (1957) 1119. 

[4] ANDERSON M.H. et al. Science 269 (1995) 198; 

BRADLEY C.C. et al, Phys. Rev. Lett. 75 (1995) 1687; 
DAVIS K.B. et al, Phys. Rev. Lett. 75 (1995) 3969; 
FRIED D.G. et al, Phys. Rev. Lett. 81 (1998) 3811; 
ANGLIN J.R. and KETTERLE W., Nature 416 (2002) 211. 

[5] KRAUTH W., Phys. Rev. Lett. 77 (1996) 3695; 

GRUTER P., CEPERLEY D., and LALOE F., Phys. Rev. Lett. 79 (1997) 3549; 

SOKOL P.E., in Bose-Einstein Condensation, ed. GRIFFIN A., SNOKE D.W., STRINGARI S., (Cam- 
bridge, Cambridge) 1995. 

[6] LEE T.D. and YANG C.N., Phys. Rev. 112 (1958) 1419. 



8 



[7] FETTER A.L. and WALECKA J.D., Quantum Theory of Many-Particle System, (McGraw-Hill, New 
York) 1971. 



HUANG K., Statistical Mechanics, (John- Wiley & Sons, New York) 1963. 
KAHN B. and UHLENBECK G.E., Physica 5 (1938) 399. 

PATHRIA R.K., Statistical Mechanics, 2nd Ed., (Butterworth-Heinemann, Oxford) 1996. 
DAI W.-S. and XIE M., Phys. Lett. A 311 (2003) 340. 
DAI W.-S. and XIE M., Phys. Rev. E 70 (2004) 016103. 
KAC M., Am. Math. Monthly 73 (1966) 1. 
DAI W.-S. and XIE M., in preparation. 

KLAIN D.A. and ROTA G.-C., Introduction to Geometric Probability, (Cambridge, Cambridge) 1997. 

BERNSTEIN R.B., ed., Atom-Molcculc CoUision Theory, (Plenum, New York) 1979. 

BAYM G., BLAIZOT J.-P., HOLZMANN M., LALOE F., and VAUTHERIN D., Eur. Phys. J. B 24 (2001) 
107. 

WU T.T., Phys. Rev. 115 (1959) 1390. 

We would like to thank the anonymous referee of our paper Europhys. Lett., 72 (6), pp. 887-893 (2005). 
Before he or she helpfully points out that the differential equation of N is solvable, we do not know this 
due to the complex form of this equation. 



9 



